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Abstract. We investigate the effect of the quark-gluon to hadron transition on the evolution of 
cosmological perturbations. If the phase transition is first order, the sound speed vanishes during 
the transition, and density perturbations fall freely. The primordial Harrison-Zel'dovich density 
fluctuations for scales below the Hubble radius at the transition develop peaks, which grow linearly 
with the wavenumber, both for the hadron-photon-lepton fluid and for cold dark matter. The large 
peaks in the spectrum produce cold dark matter clumps of 10~ 8 to 1O _11 M0. 

QCD makes a transition from a quark-gluon plasma at high temperatures to a hadron gas 
at low temperatures. Lattice QCD simulations give a transition temperature T* ~ 150 MeV 
and indicate a first-order phase transition for the physical values of the u,d,s-quark masses [|TJ. 
The relevance of the QCD transition for cosmology, especially for big-bang nucleosynthesis 
[0], has been discussed before, but the focus was on effects of bubble formation. In this paper 
and in || we look at matter averaged over scales much larger than the bubble separation. 
We show that for a first order phase transition the sound speed c s = {dp/dpYJ 2 drops 
to zero suddenly at the moment the transition temperature T* is reached, stays zero for 
the entire time until the phase transition is completed, and afterwards suddenly rises back 
to c s ~ c/ v^3- In contrast the pressure stays positive and varies continuously, although 
it goes below the radiation fluid value p/p = 1/3. Since c s is zero during the transition, 
there are no pressure perturbations, no pressure gradients, no restoring forces. Pre-existing 
cosmological perturbations, generated by inflation with a Harrison-Zel'dovich spectrum, go 
into free fall for about a Hubble time. The superhorizon modes (at the time of the transition) 
remain unaffected, the subhorizon modes develop peaks in Sp/p which grow linearly with 
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wavenumber, ~ k/k+, where k^ hjs ~ Hubble rate H at the end of the QCD transition. 

The sound speed, c s = (dp/dpYJ 2 , must be zero during a first-order phase transition of 
a fluid with negligible chemical potential, since the fluid must obey 

> + P = T%, (!) 

according to the second law of thermodynamics. Because the energy density p is discontinu- 
ous in temperature at T* for a first-order phase transition, the pressure p must be continuous 
with a discontinuous slope. As the universe expands at fixed temperature T* during the phase 
transition, p as a function of time slowly decreases from p+(T*) to p_(T*), p stays constant 
at p(T^), and therefore c s is zero during the whole time of the phase transition. 

The interaction rates in the QCD-photon-lepton fluid are much larger than the Hubble 
rate, T/H ^> 1, therefore we are very close to thermal and chemical equilibrium, the QCD 
transition is very close to a reversible thermodynamic transformation. Estimates show that 
supercooling, hence entropy production, is negligible, (T* — T supcrcoo i ing )/T^ ~ 10~ 3 0. Bub- 
ble formation is unimportant for our analysis, estimates give a bubble separation ~ 1 
cm ||, while the Hubble radius at the QCD transition is Rh ~ 10 km. We shall analyze 
perturbations with A 3> 

The bag model gives a parameterization and a reasonable fit to the lattice QCD data ||. 
In the bag model it is assumed that for T > T* the quark-gluon plasma (QGP) obeys 

PqgAT) = P%gf(T) - B , (2) 

where Pqop(T) = (7t 2 /90)(7q GP T 4 , g* is the effective number of relativistic helicity states, 
and B is the bag constant. We include u,d-quarks and gluons in the quark-gluon plasma, 
and for T < T* we have a hadron gas (HG) of pions, which we treat as massless and ideal, p 
follows from Eq. (||D via the second law, Eq. (JJ), and s from s = dp/dT. The bag constant is 
determined by the critical temperature T* via pqgp(T^) = pn G (T*). The photon- lepton fluid 
contains 7, e, /1, and 3 neutrinos. The growth of the scale factor during the QCD transition, 
a + /a_ ~ 1.4, follows from the conservation of entropy in a comoving volume. 

The evolution of linear cosmological perturbations through the QCD transition is ana- 
lyzed in the longitudinal sector (density perturbations) for perfect fluids. We choose a slicing 
S of space-time with unperturbed mean extrinsic curvature, S[trKij(T,)} = 0. The adapted 
gauge is the uniform expansion (Hubble) gauge 0. As fundamental evolution equations for 
each fluid we have V^T^ = 0, i.e. the continuity equation and (in the longitudinal sector) 
the 3-divergence of the Euler equation of general relativity, 

d t e = -3#(e + vr) - Aip-3H(p + p)a (3) 
d t tp = -3Hi; - 7T - {p + p)a , (4) 

where e = 5p, n = 6p, p = po, p = po, Vip = momentum density, a = lapse function. The 
system of dynamical equations is closed by Einstein's i^g-equation, the general relativistic 
Poisson equation, 

(A + 3H)a = 47rG(e + 3tt) , (5) 
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Figure 1: The modifications of the CDM density contrast A CDM = |<5cdm| at T+/1Q and 
of the radiation fluid amplitude A R ad = (^rad + ^rad) 1 ^ 2 due to the QCD transition. 
Both quantities are normalized to the pure Harrison-Zel'dovich radiation amplitude. On the 
horizontal axis the wavenumber k is represented by the CDM mass contained in a sphere of 
radius n/k. 

together with the equation of state. These three equations are the Jeans equations extended 
to general relativity in the longitudinal sector. 

Our numerical results of evolving a mode k of a cosmological perturbation through the 
QCD transition are given in Fig. [1|. We work with the dimensionless variables 5 = e/p 
(density contrast) and ip = k phys ip/p (~ peculiar velocity). We show the enhancement of the 
amplitude A R ad = (^rad + 3^rad) 1//2 °f the acoustic oscillations of the radiation fluid (QCD, 
photons, leptons) after the transition compared to the amplitude without transition. For 
cold dark matter (CDM) we show the amplitude Acbm = I^cdmI at T*/10 compared to A R ad 
without transition. In both cases we obtain peaks over the Harrison-Zel'dovich spectrum 
of adiabatic density fluctuations. Only those modes are affected which are subhorizon at 
the end of the transition, k>k±, where k phys (t + ) ~ H(t + ). Our peaks grow linearly in k 
for k 3> k+. The modes k are labeled by the CDM mass contained in a sphere of radius 
A/2 = ir/k. The value k* corresponds to Mq DM ~ 1O _8 M . The radiation energy inside 
A*/2 is ~ 1M , but it gets redshifted as Mrad(«) ~ (a e quaiit y /a)^fcDM- The perturbations 
in the radiation fluid will be wiped out by collisional damping from neutrinos on scales much 
smaller than the horizon scale at neutrino decoupling (1 MeV). 

For cold dark matter we consider the lightest supersymmetric particle (LSP), the neu- 
tralino Xi m t ne minimal supersymmetric standard model with universal supersymmetry 
breaking at the grand unified scale. LEP 1.5 combined with the gluino mass limit from 
Fermilab gives a minimal mass M£gp = 27 GeV 0. With a freeze-out temperature of 
Tf ~ Mlsp/20 free-streaming wipes out CDM structure for k/k± > 10. 
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The origin and magnitude of these impressive peaks for > ^ is easily understood. 
The radiation fluid in each subhorizon mode makes standing acoustic oscillations before and 
after the QCD transition with gravity negligible and with equal amplitudes of S and \[?>i\). 
During the transition the sound speed is zero, there are no restoring forces from pressure 
gradients, the radiation fluid goes into free fall. If the transition time is short, At < if -1 , 
gravity is again negligible for the radiation fluid during this free fall. This is inertial motion 
in the sense of Newton. The peculiar velocity is constant, and the density contrast grows 
linearly in time with a slope k. Therefore the height of the peaks is (v4+/v4_) pca ks ~ k/k±. 
Hence modes with k phys / H>10 4 go nonlinear by the end of the QCD transition. The explicit 
solution along these lines is presented in ||. CDM falls into the gravity wells generated by 
the radiation fluid. 

The peaks above the Harrison-Zel'dovich spectrum lead to CDM clumps with 10~ 8 > 
Mqdm/Mq > 10 -11 , which go nonlinear sometime after equality and virialize by violent 
gravitational relaxation. The mass range of these CDM clumps lies just beyond the smallest 
mass (5 x 1O _8 M ) accessible to present microlensing observations ||, thus future lensing 
observations might discover such clumps. 
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